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ABSTRACT 

We investigate departures from local thermodynamic equilibrium (NLTE) in the line forma- 
tion of neutral and singly ionised iron lines and their impact on spectroscopic stellar parame- 
ters. The calculations are performed for an extensive grid of ID MARCS models of metal-rich 
and metal-poor late-type dwarfs and giants. We find that iron abundances derived from Fel 
lines are increasingly underestimated in hotter, lower surface-gravity, and more metal-poor 
stars, in a simple and well-defined pattern, while LTE is usually a realistic approximation for 
Fe II lines. For the vast majority of dwarfs and giants, the perturbed ionisation balance of Fe I 
and Fe II is the main relevant NLTE effect to consider in the determination of spectroscopic 
parameters, while for extremely metal-poor stars and hot giant stars significant impact is seen 
also on the excitation balance and on the microturbulence determination from Fe I lines. 
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(1D) geometry. The most critical approximation is the mixing- 
length description of the convective energy transport. The advent 
of 3D, radiation-hydrodynamical simulations has shed light on the 
systematic uncertainties introduced by these simplifying assump- 
tions. In particular, metal-poor stars have been demonstrated to 
hav e much cooler line-forming layers than previously thought (see 
e.g. lAsplundll2005l and references therein). 

In Paper I in this series dBergemann et alj|2012h we demon- 
strated that a modelling technique allowing for departures from 
LTE can be used to accurately predict iron abundances and spec- 
troscopic stellar parameters for a set of benchmark late-type stars. 
We specifically illustrated how traditional ID hydrostatic models 
successfully meet the ionisation and excitation equilibria for solar 
metallicity stars, but evidently fail to achieve excitation equilibrium 
at realistic temperatures in metal-poor stars. In particular, low exci- 
tation lines are sensitive to the atmospheric structure, and require a 
combination of 3D and NLTE analysis. We further illustrated how 
the ionisation balance of Fe I and Fe II lines can be exploited to 
infer realistic spectroscopic surface gravities and metallicities also 
with traditional ID models, once the effective temperature has been 
determined by other means. Alternatively, the ionisation balance 
can be used to constrain the effective temperature if the gravity is 
known by other means. 

Due to the increasing numerical complexity, compared with 
the LTE case, NLTE investigations have previously been limited to 
individual stars and usually only a handful of spectral lines. For the 
first time, we present an extensive grid of calculations, with indi- 
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1 INTRODUCTION 

Measuring the iron content of late-type stars is one of the primary 
challenges for the field of Galactic archaeology. Due to its large 
opacity contribution in cool stellar atmospheres, Fe has come to 
serve as a fundamental reference point for all chemical analysis and 
its interpretations. Furthermore, the visibility of a wealth of spec- 
tral lines with a range of atomic properties and line strengths en- 
ables the determination of spectroscopic stellar parameters through 
the excitation and ionisation equilibria. With the knowledge of the 
star's iron content, effective temperature, and surface gravity, we 
can then acquire further information about the evolution of the star 
itself, i.e. determine its mass and age, and the evolution of the stel- 
lar population it resides in. 

However, the low fraction of neutral iron in stellar atmo- 
spheres makes its line formation sensitive to departures from lo- 
cal thermodynamic equilibrium (LTE). Traditional LTE analyses 
of Fel lines therefore tend to underestima te the true Fe abundance , 
as is well described in the literature ( e.g. Thevenin & Idiardl 1 9991 : 



iGehren et alj200ll : ICollet et alj|2005l : iMashonkina et al.ll201ll)~ 

The success of spectroscopic methods clearly depends on the 
realism of the atmospheric structure. Even if high precision may be 
obtained with a very simplistic model, accuracy may or may not. 
Traditional methods rely on model atmospheres calculated under 
the assumptions of hydrostatic equilibrium and a one-dimensional 
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Table 1. Dimensions of model atmosphere grid 
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Figure 1. Example LTE and NLTE curves-of-growth for a UV Fe I tran- 
sition. The model parameters are 7" e (j = 6500 K, log(g) = 4.0, and 
^, = 2km/s. The NLTE abundance correction at a given equivalent width, 
AA(Fe)i, is defined as the difference between the two curves. 



vidual NLTE abundance corrections given for thousands of lines in 
the ultra-violet, optical and near-infrared part of the spectrum. The 
calculations span from ultra metal-poor to super-solar metallicities 
and include dwarfs, subgiants, and red giant stars. The calculations 
can easily be extended to an analogous grid of spatially and tem- 
porally averaged 3D models (her eafter <3D> ), onc e such becomes 
available from e.g. the Stagger dCollet et alj|201 lh or C05BOLD 
collaborations (iLudwig et alj|2009l) . 



2 METHODS 

As is customary for late-type stars, we solve the restricted NLTE 
problem given the assumption of iron as a trace element. Some jus- 
tification for this methodology is presented in Sec t.| 3.3.3l The ID, 
NLTE code MULTI, version 2.3, jCarlssonlll98rj,|l992h was used 
to simultaneously solve the equations of radiative transfer and sta- 
tistical equilibrium for a model atom of neutral and singly ionised 
iron. 
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tion were computed using semi - empirical reci pes l lvan RegemorteJ 
1 1 9621 ; ISeator]! 1 9621 ; I AllerJ 1 9731 ; iTakedall 1 994h . 

The efficiency of excitation and ionisation (and charge ex- 
change reactions) induced by collisions with neutral hydrogen is 
commonly the largest source of uncertainty in statistica l equilib- 
rium calculations for late-type stars (see discussion in e.g. lAspl und 
120054) . In the absence of better alternat ives, we adopt the Drawin 
formulae dSteenbock & Holwegen i ri9840 for these processes, and a 
free scaling parameter (Sh) that sets the absolute efficiency of the 
collisions. As explained in Paper I, the scaling factor can be con- 
strained when confronted with observations of stars with surface 
gravities and effective temperatures accurately known through in- 
dependent methods. We discuss the impact of our choice of scaling 
factor, S H = 1.0, in Sect. |3~3~2l 

When the atomic level populations had reached convergence, 
high-resolution NLTE line profiles of 3239 lines of Fel and 107 
lines of Fell were computed. The lines were selected based on 
having a non-negligible line strength at solar metallicites and all 
fall in the wavelength range 3200 - 9300A . We impose no fur- 
ther constraints based on the accuracy of atomic data or contri- 
bution of blending species in stellar spectra. Hence, lines may or 
may not be suitable for abundance a nalysis in a giv en star. The os- 
cillator strengths are adopted from dKuruczl 120071) . For 1479 Fel 
lines and 95 Fell lines, ABO theory and the quantum mechan- 
ical treatment of pressure broadening by H collisions is imple- 
mented for broadening by neutral hydr ogen dBarklem etaflbOOCt 
iBarklem & A spe und-Johansson 2005j). In the absence of these 
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data, we resort to 
of 1.5. Stark broadening is neglected. 



2.1 Model atom 

The model atom used to establish the statistical equilibrium was 
presented in detail in Paper I and we give only a summary here. 
The atom consists of 296 levels of Fe I, 112 levels of Fe II and the 
Fe III ground state, with experimentally determined and/or theoret- 
ically predicted level energies. Levels close in energy are collapsed 
into superlevels to save computing time. Radiati ve bound-bound 
transi tion probabilities are collected from NISlQ l lRalchenko et alj 
l2012h and R. L. Kuruc^l data bases a nd pho to-ionisation cross- 
sections were provided by M. Bautista d 19971 and priv. communi- 
cation). Rate coefficients for electron impact excitation and ionisa- 

1 http://physics.nist.gov/asd 

2 http : //kurucz . harvard . edu/| 



2.2 Model atmospheres 

We use a grid of ID, h ydrostatic model atmos pheres computed 
with the MARCS code dGustafsson et alj|2008l). The models al l 
have standard composition^!, i.e. scaled solar dGrevesse et all2007tl 
with gradual alpha-enhancement up to 0.4 dex at [Fe/H] < -1.0. 
The grid dimensions are specified in TableQ] Models that are miss- 
ing from the original grid are obtained by interpolation or, in some 
cases, substitutiorQ. 



3 See lmar c s.astro.uu.sel 

4 When needed and justified we do not enforce consistency between the 
microturbulence used to compute the model atmosphere and that used to 
perform the NLTE calculations, i.e. missing models are replaced with ones 
of higher/lower microturbulence velocities. 
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Figure 2. The lines mark the typical sizes of NLTE effects for high-excitation (£ C xc > 2.5 eV), unsaturated (W,\ < 50 mA ) Fe I lines and their dependence on 
stellar parameters. All models shown have f t = 2.0km/s. 



2.3 NLTE abundances 



Here we base the comparison between LTE and NLTE on equiva- 
lent widths (Wa), rather than synthetic line profiles. Note however 
that differences in lines profiles become important for saturated 
lines observed at high resolution. For each grid point, we compute 
the LTE and NLTE equivalent widths for the iron abundance that is 
identical to the one adopted for the atmosphere. The line strengths 
increase smoothly with increasing metallicity in a manner similar 
to a traditional curve-of-growth. Fig.[T]illustrates the curves derived 
for a spectral line at a certain effective temperature, surface gravity, 
and microturbulence. The NLTE abundance correction is defined as 
the difference between the NLTE and LTE curves at a given equiv- 
alent width. For practical reasons we introduce a cut towards weak 



lines at 1 mA . Lines weaker than this are not discussed. Our equiv- 
alent width integrations are numerically reliable up to ~1000 mA . 

Inspection of the NLTE abundance corrections determined for 
a given set of stellar parameters reveals a strong similarity between 
lines of similar excitation potential and line strength. In the follow- 
ing discussion we will use the concept of a mean NLTE abundance 
correction, restricting ourselves to all lines below a certain limit in 
equivalent width (here 50 mA ) and above a certain limit in excita- 
tion potential (here 2.5 eV) of the lower level. The first restriction 
is made because abundances of strong lines are through common 
practice forced into agreement with weaker lines by fine-tuning of 
the microturbulent velocity. The absolute abundances and NLTE 
corrections derived from microturbulence-sensitive lines are thus of 
lesser importance. The second restriction to high-excitation lines is 
motivated by their lower sensitivity to 3D effects compared to low- 
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excitation lines, as discussed in Paper I. However, in ID, our atomic 
model predicts only minor differential NLTE effects for lines of 
high and low excitation potential. Figure[2]illustrates how the mean 
NLTE abundance correction of un-saturated, high-excitation Fel 
lines (AA(Fe)i) varies with stellar parameters. 



3 DISCUSSION 

3.1 Dependence on stellar parameters 

As Fig.[2]shows, the departures from LTE have a smooth and simple 
dependence on stellar parameters. In summary, NLTE abundance 
corrections increase with decreasing metallicity, increasing effec- 
tive temperature, and decreasing surface gravity. Here we do not 
aim to disentangle all the contributing factors in order to explain 
the detailed behaviour quantitatively, but rather we highlight some 
apparent circumstances that enable a qualitative understanding of 
the trends. 

We focus here on the perturbed ionisation balance, established 
by over-ionisation in Fe I bound-free transitions compared to LTE. 
We take the term overionisation to imply that - njPji > 0, 
where n, is the level population of a b ound state of F e I, and is the 
ground state of Fe II (notations follow lRuttenl2003l) . Py = i?y + Cy, 
where R t j is the radiative rate from level i to j and Cy the cor- 
responding collisional rate. Hence, over-ionisation means that the 
ionisation-recombination balance for a transition is shifted towards 
more efficient ionisation, which tends to de-populate the lower level 
with respect to LTE. Analogously, over-excitation is adopted for a 
transition whose excitation-de-excitation balance is shifted towards 
more efficient excitation. Fig.[3] demonstrates how the total photo- 
ionisation and photo-recombination rates per unit volume (n,/? y and 
njRji) vary with stellar parameters and the threshold wavelength of 
the transition. Also shown are the total collisional rates, which in 
LTE are related by n,C,j = njCp. 

The photo-ionisation rate, /Jy, is governed by the photo- 
ionisation cross-sections and the radiation field. In practice, the 
second dominates the energy dependence of /?y, which exhibits a 
sharp increase from UV to optical wavelengths and thereafter de- 
clines mildly towards infra-red wavelengths |j. The Fel level popu- 
lation is maximal for the ground state and decreases exponentially 
with excitation energy. The product ra,/fy thus behaves as shown in 
Fig-El peaking in the UV regime. The recombination rates per unit 
volume (njRjj) follow a very similar behaviour but deviate from 
the photo-ionisation rates as the radiation field deviates from the 
Planck function. 

Collisional rates, Cy, on the other hand, depend on the bound- 
free collisional cross-sections and the availability of impact elec- 
trons and atoms. The impact species have kinetic energy spectra 
that peak at energies corresponding to infra-red wavelengths □ . In 
addition, the bound-free collisional cross-sections for electrons and 
hydrogen atoms increase with decreasing energy gap. These two 
dependencies outweigh the drop in level population, and the total 
rate uniformly increases towards longer wavelengths. The relative 
behaviour of the plotted curves is tightly connected to the sizes of 

5 Low flux levels in the far UV in late-type stars suppresses the photo- 
ionisation rates for low-excited levels. Rates for transitions that bridge 
smaller energy gaps grow increasingly larger as they feel the effect of the 
stronger radiation field in the optical regime. 

6 At typical temperatures, kT ~ 0.5 eV, which corresponds to a kinetic 
energy peak at 2/im. 



the departures from LTE, through the equations of statistical equi- 
librium. 

Of course, all levels are also coupled by bound-bound radia- 
tive rates that can drive the level populations from LTE. Lines, or 
parts of lines, which are formed close to continuum-forming lay- 
ers (t 50 o ~ 1) in the atmosphere react analogously to bound-free 
transitions to the non-Planckian radiation field. Excitation is thus 
favoured above de-excitation in the UV and optical in the inner at- 
mosphere, more specifically in layers which are optically thick to 
the radiation in strong line cores. The accumulated effect of the 
imbalance in the bound-bound rates is here to enhance the under- 
population of relevant levels of neutral iron. This can be understood 
because the levels most affected by over-excitation are the lowest 
excited, while those most affected by over-ionisation are moder- 
ately excited (2-4 eV, see Fig.[3](- The level populations of Fel 
are thus re-distributed by line transitions in a way that enhances 
under-population of relevant Fel levels (< 5 eV). We note that in 
shallow atmospheric layers, where the cores of strong lines become 
optically thin, different NLTE effects set in that instead favour de- 
excitation. 

We now proceed to interpret the NLTE dependencies on stel- 
lar parameters. The over-ionisation is driven by the surplus of UV 
photons with respect to the Planck function in the atmospheric radi- 
ation field. In turn, the size of the J v - B v -excess in the UV depends 
on the steepness of the atmospheric temperature gradient and the 
efficiency of metal line blocking in this part of the spectrum. We can 
therefore easily interpret the trend of increasing NLTE effects with 
decreasing metallicity. It can be seen in Fig.[3]that the gap between 
the photoionisation and photo-recombination curves grows larger 
at lower metallicity, while the relative contribution of radiative and 
collisional processes is approximately constant (the absolute rates 
are obviously lowered in accordance with the smaller amount of Fe 
atoms). 

The clear trend of increasing NLTE effects at lower surface 
gravity can be understood from Fig.[3]as due to the smaller relative 
strengths of collisional rates, scaling with the number density of 
iron atoms and the ionising impact species, with respect to radiative 
rates, that only scale with the number density of iron. 

The temperature gradient steepen sharply at higher effec- 
tive temperatures, while the flux maximum moves towards shorter 
wavelengths. In Fig.[3] we see that higher effective tempera- 
tures again widen the gap between photoionisation and photo- 
recombination curves. In addition, the relative efficiency of colli- 
sions with respect to radiative transitions in lessened, since the in- 
creased amount of photons outnumber the increased amount and 
higher speed of the impact species. All in all, this explains why 
higher effective temperatures give rise to larger NLTE effects. 

Finally, one might ask why only the strengths of Fel lines 
are affected by NLTE and not Fell lines, for which similar ar- 
guments as outlined above hold at least qualitatively. Bound-free 
transitions in Fe II are much less influential on the statistical equi- 
librium because of the large ionisation potential, hence low num- 
ber density of Fe III. Furthermore, the three stellar-parameter de- 
pendencies (lower metallicity, lower surface gravity, and higher ef- 
fective temperature) that give rise to increasing influence from the 
UV radiation field also act as to raise the relative fraction of singly 
ionised iron, which usually is the dominant species. A perturba- 
tion from the LTE populations clearly has a larger relative influ- 
ence on a minority species. The well populated levels from which 
the visible Fell transitions originate (< 4eV) are generally close 
to thermalised throughout the whole atmosphere, while departures 




Figure 3. The panels show the total number of Fe I-Fe II ionisation and recombination events per unit time and volume for selected models at line forming 
layers (log(tau5oo = _ 2)). The abscissa marks the threshold wavelength, n\ is the LTE level population of the lower level, and nj the LTE population of 
the upper level of the corresponding transition, i.e. the ground state of Fe H. The curves have been constructed by applying a smoothing factor of ten to the 
individual values, so that each point on the curve represents the average of ten adjacent rates. 



from LTE are apparent only above very high excitation energies 
(> 8eV). 



3.2 Determining stellar parameters 

In the following subsection we will discuss the differences between 
LTE and NLTE stellar parameters as inferred from Fe I and Fe II 
lines. 



3.2.1 Metallicity 

How NLTE affects the metallicity determination can be straight- 
forwardly estimated from Fig.|2] Spectroscopic LTE analyses un- 
derestimate the iron abundance derived from Fe I lines by anything 
between 0.0 and 0.5 dex, the exact value being strongly dependent 
on T eff , logg and [Fe/H], Abundances derived from Fell lines are 
on average not affected at all in models with [Fe/H] > -3. Even 
at [Fe/H] = -4, the corrections for the few lines that are still 
detectable are smaller than 0.02 dex. We conclude that LTE is a 
safe assumption for Fe II lin es in our parameter space, in line with 
severa l previous studies, e .g. Thevenin & Idiart ( 1999): lKorn et alj 
fe003h:lcollet et alj j2005t) : lMashonkina et alj d201 lh . but see also 
ICrametalJ l ll98fj|) 



3.2.2 Surface gravity 

It is common to vary the model surface gravity until agreement 
is achieved between abundances derived from neutral and singly 
ionised iron lines, so called ionisation balance. Since Fe I lines are 
much less sensitive to surface gravity variations compared to Fe II 
lines, the metallicity is then practically based on the former. Since 
Fe I lines are subject to significant NLTE effects, the LTE ionisation 
balance is not always realistic and the surface gravity as well as the 
metallicity will therefore be underestimated with this method. We 
have assessed the size of this effect by interpolating LTE curves- 
of-growth with varying surface gravity onto the NLTE equivalent 
widths found at each grid point. Since we are only interested in the 
difference between LTE and NLTE, the comparison can in princi- 
ple be inverted and the LTE values chosen as the reference values. 
This would give very similar results, but since the NLTE equivalent 
widths are presumably more realistic it is more appropriate to use 
them as reference. For this comparison we again limit ourselves to 
un-saturated, high-excitation Fel and Fell lines. 

The results, i.e. Alogg = logg NLTE - logg LTE , are shown in 
Fig.|4] As a rule-of-thumb, the surface gravity is underestimated by 
a factor x times the size of the NLTE correction for a given model, 
in logarithmic units. For dwarfs, x ss 2.5 while for giants the sen- 
sitivity increases slightly, x ~ 3. Thus for a metal-poor turn-off 
star at T eS = 6500 K, log(g) = 4.0, [Fe/H] = -2, that has a typi- 
cal NLTE abundance correction of 0.1 dex, LTE analysis underes- 
timates the surface gravity by 0.25 dex. An important consequence 
of this behaviour is that metal-poor dwarfs easily can be misclas- 
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Figure 4. The lines mark the estimated NLTE effect on surface gravities derived from ionisation balance of high-excitation, unsaturated Fe I and Fe II lines, 
given that other stellar parameters are fixed by independent means. All models shown have f t = 2.0km/s. 



sified as subgiants. Similarly, a red giant star with r eff = 4500 K, 
log(g) = 1.5, [Fe/H] = -5 that has a NLTE abundance correction 
of 0.2 dex would have its surface gravity underestimated by 0.6 dex 
in LTE. 

3.2.3 Effective temperature 

Here, we consider two different methods to infer the effective tem- 
perature, while the other parameters are kept fixed. First, given a 
reliable measurement of surface gravity Q one may constrain ef- 
fective temperature through ionisation balance. In this case, Fel 



based abundance are effectively brought into agreement with the 
less temperature sensitive Fell lines. The metallicity determina- 
tion itself thus does not suffer heavy bias in the LTE case, but 
the optimised effective temperature will be overestimated, by the 
amount shown in Fig|5] The calculations were based on weak, high- 
excitation lines, as before. 

Secondly, the effective temperature can be inferred from the 
excitation balance of Fe I lines, flattening abundance trends with 
the excitation potential of the lower level of the transition (Fell 
lines typically do not span a sufficient range in excitation potential 
in spectra of late-type stars). The method exploits the larger temper- 



7 This can be provided by an accurate parallax measurement or from stellar 
evolution calculations. Both methods require a reasonable estimate on ef- 



fective temperature, in which case the stellar parameter determination may 
need iteration. 




Figure 5. The lines mark the estimated NLTE effect on effective temperatures derived from ionisation balance of high-excitation, unsaturated Fe I and Fe II 
lines, given that other stellar parameters are fixed by independent means. All models shown have f t = 2.0km/s. 



ature sensitivity of low excitation lines compared to high excitation 
lines, which simultaneously makes them more sensitive to depar- 
tures from LTE. For ID models, NLTE abundance corrections are 
nevertheless only mildly differential with excitation potential, but 
they become very significant for metal-poor <3D> models, as il- 
lustrated in Paper I. 

We have estimated the difference between ID LTE and NLTE 
effective temperatures derived through excitation equilibrium, by 
again interpolating the LTE curves-of-growth onto NLTE results 
for all grid points, and enforced a flat trend with excitation poten- 
tial. Again, we restrict ourselves to unsaturated lines, but the only 
restriction made in excitation potential is that at least one transi- 
tion must originate from above 3.5 eV. Fig. [6] shows the typical 
difference in effective temperature derived in the two cases, i.e. 
Ar c ff_ cxc = T e (f_NLTE - r e (f LTE- Because NLTE abundance corrections 



are generally more positive for low excitation lines, the effective 
temperatures derived through LTE analysis are usually higher. Not 
surprisingly, the effect tends to increase in tandem with the NLTE 
effects, but the shape of curves are somewhat different. Down to 
metallicities as low as [Fe/H] = -3 the curves are rather flat and 
all types of dwarfs and the bulk of red giant branch stars are af- 
fected by less than 50 K. Excitation temperatures derived with ID 
models are thus not heavily influenced by NLTE effects except for 
horisontal branch stars and extremely metal-poor stars. However, 
we caution that even ID NLTE results are not reliable in situations 
where 3D effects are pronounced, in particular if departures from 
LTE are minor. 
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Figure 6. The lines mark the estimated NLTE effect on effective temperatures derived from excitation balance of un-saturated Fel lines, given that other stellar 
parameters are fixed by independent means. All models shown have f t = 2.0km/s. 



3.2.4 Microturbulence 

The NLTE abundance corrections show a dependence on line 
strength beyond saturation (W A > 50 mA), which will impact on 
the determination of microturbulence from Fe I lines. While the ef- 
fect may be quite significant for individual lines, the typical change 
is very small. As a general rule, saturated (50 < W A < 100 mA) 
lines have an abundance correction that is larger than for un- 
saturated lines, the difference being roughly factor of 5 less than the 
NLTE abundance correction itself in logarithmic units. Ultimately, 
this leads to lower microturbulent velocities in LTE compared to 
NLTE. 

The size of this effect is illustrated in Fig.[7] These results 
have been obtained analogously to the other parameters by calcu- 
lating the microturbulence value that flattens the LTE abundance 



dependence on reduced equivalent width, using the NLTE equiv- 
alent widths as reference in each grid point. The calculations are 
restricted to high-excitation lines, to minimise a systematic bias as 
explained above. We also introduce an upper cut in line strength 
at log(W,{//0 = -4.6 to exclude the damping part of the curve- 
of-growth. Similarly to the excitation balance, the effect is essen- 
tially negligible for dwarfs and giant stars above [Fe/H] > -3. Ex- 
tremely metal-poor stars are more affected, but generally they have 
few lines that are sensitive to this parameter. 
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Figure 7. The lines mark the estimated NLTE effect on microturbulence velocities derived from high-excitation Fe I lines, given that other stellar parameters 
are fixed by independent means. The reference value in LTE is 1 km s _I 



3.3 Modelling uncertainties 

3.3.1 Model atmospheres 

Accurate absolute abundances in LTE and NLTE require realistic 
model atmospheres, hence spectrum synthesis based on 3D hydro- 
dynamical models are expected to give superior results. Our first 
attempts to investigate this problem in LTE and NLTE using the 
averages of the 3D simulations of stellar convection were illus- 
trated in Paper I for selected well-observed stars. We found that the 
average metallicities derived from <3D> models are not dramati- 
cally different from ID models, in particular in NLTE (< 0.04 dex). 
Most sensitive to the model structure are abundances derived from 
saturated, low excitation lines in metal-poor stars. The greater sim- 
ilarity between NLTE abundances of ID and <3D> models can be 
understood simply via the close resemblance of the mean radiation 



field throughout the atmospheres, as it is largely set by conditions at 
continuum optical depth 7"5oonm ^ L contrary to the Planck function 
that behaves very differently in shallow atmospheric layers. Paper 
I also illustrates how stellar parameters inferred through ionisation 
equilibrium of Fe I and Fe II lines agree well with NLTE modelling 
in ID and <3D>. We emphasize, however, that the excitation equi- 
librium established with <3D> models is generally superior to that 
of ID models. Finally, it is worth to note that the uncovered depen- 
dencies of the departures from LTE on the stellar parameters are 
likely to hold true for any type of model atmosphere grid. 

3.3.2 Atomic data 

Our analysis focuses on systematic differences between LTE and 
NLTE, which are not affected by the choice of oscillator strengths 
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and broadening data adopted for spectrum synthesis. According 
to our assessment, atomic data for level energies, radiative tran- 
sition probabilities, and photo-ionisation cross-sections are suffi- 
ciently advanced and complete not to introduce large uncertainties 
in the modelling. Also, according to our tests the remaining uncer- 
tainties in the electron impact collision rates do not to introduce 
large uncertainties in the NLTE modelling. As mentioned in the 
introduction the most critical parameter is the efficiency of colli- 
sions with neutral hydrogen, which we constrained empirically in 
Paper I. Lowering the scaling factor Sh by one order of magni- 
tude, from 1.0 to 0.1, increases the NLTE effects by typically a 
factor 2-4, depending on the stellar parameters. The qualitative 
behaviour remains the same, but the influence of the radiative-rate 
imbalance that causes over-ionisation increases. Paper I demon- 
strates that such large NLTE effects are unrealistic for metal-poor 
dwarfs and subgiant stars, but are potentially more successful for 
the metal-poor giant HD122563. Here, we have chosen to calculate 
the NLTE grid for Sh = 1.0, the most supported value, only. 

3.3.3 The trace element assumption 

As mentioned in the introduction we solve the restricted NLTE 
problem, which neglects all potential feedback that the perturbed Fe 
level populations have on the atmospheric structure. This simplifi- 
cation is necessary to make the computations tractable and the con- 
sequences have not been rigorously investigated. The most note- 
worthy implicit assumptions are that the demonstrated changes in 
excitation and ionisation balance have no significant implications 
for the free electron density and the UV bound-bound and bound- 
free opacity. Given the known size of the LTE departures over the 
late-type star grid, we now proceed to make simple tests of the ex- 
tent of such a feedback. 

The main opacity change would be caused by the lower 
amount of neutral iron compared to LTE, which in practice can be 
mimicked with an atmosphere of lower metallicity. Furthermore, 
NLTE populations of iron would provide more free electrons in 
the atmosphere. The increase in electron number density can be 
straightforwardly estimated by the decrease in neutral iron, since 
each additional iron ion corresponds to one additional electron. It 
turns out that this effect is negligible over the whole stellar grid, 
since the free electron density is typically much larger than that of 
neutral iron (hydrogen and metals with lower ionisation potential 
are the main electron donors). In the most extreme case, the elec- 
tron density of a cool, metal-rich giant would change by < 1%, 
which is insignificant in this context. We therefore disregard the 
second aspect and focus our test on the opacity change. 

We selected a few models from the grid and performed new 
NLTE calculations on models with the same stellar parameters, 
but with a metallicity lower by the mean NLTE correction. The 
lower metallicity was also adopted for the calculation of continuous 
and line background opacity. Thereafter we compare the LTE and 
NLTE equivalent widths for a given abundance, obtained from the 
two different model atmospheres. We perform the test on a metal- 
rich and metal-poor dwarf and giant with the same parameters as 
shown in Fig.|3] 

The solar metallicity models have a higher sensitivity to the 
metal content and are more affected by the mentioned inconsis- 
tency, even though the predicted sizes of the NLTE corrections are 
much smaller than for metal-poor stars. In particular, continuous 
opacities in the UV decrease when the model metallicity is low- 
ered. We thus expect equivalent widths of Fe lines to increase due 
to lowered continuous opacities, while the NLTE effects are ex- 



pected to increase and weaken the Fe I lines with respect to LTE. 
In practice, the NLTE equivalent widths inferred from the model 
with r eff = 6500 K and log(g) = 4.0 increase by up to 1% when 
the model metallicity (not the abundance used for computing syn- 
thetic line profiles) is lowered by 0.03 dex, corresponding to the 
mean NLTE abundance correction. This is certainly a second-order 
effect, approximately ten times smaller than the NLTE effect itself. 

Note that our test implicitly assumes that all the ionisation bal- 
ance of all other elements are similarly affected by NLTE as Fe, 
which is not the case. Only consistent NLTE model atmosphere 
calculations will tell for sure how rea listic the trace element as- 
sumption is for Fe and other elements. IShort & Hauschildl J2005t) 
present such a model for the Sun, displaying a dramatic increase in 
UV flux in NLTE compared to LTE, and a significantly more shal- 
low atmospheric temperature structure at continuum forming re- 
gions. However, their comparison to observed absolute solar fluxes 
suggests that the opacity loss due to over-ionisation of iron-peak 
elements is over-estimated, or, alternatively, that other important 
sources of opacity are yet missing. 

3.4 Comparison with other studies 

iMashonkina et al.l d201ll) demonstrated in detail how recent devel- 
opments in atomic data calculations have improved the NLTE mod- 
elling of iron lines. In particular, R. L. Kurucz's calculations of high 
excitation energy levels beyond reach of experiments have enabled 
a realistic coupling to the next ionisation state, without having to 
enforce a n arbitrary thermalis ation above a certain energy threshold 
(see also lGehren et alj|200ll) . As described in Paper I, our NLTE 
calculations for standard stars are in g ood quantitative agreem ent 
with those of IMashonkina et alj J201 lh . In lMashonkini] d201ll) the 
author presents a small grid of calculations for A-F type dwarfs of 
solar metallicity. Their NLTE abundance corrections are of similar 
order of magnitude to ours, but differences for individual Fe I lines 
can amount to more than 0.05 dex, which is comparable to the size 
of corrections for solar metallicity stars. There could be many rea- 
sons for this, among them different choices of scaling factor of hy- 
drogen collisions. However, the dependence on surface gravity and 
effective temperature also show contradictory behaviour, which is 
surprising. While our NLTE corrections for e.g. Fel 528. 179 nm 
and 521.792 nm in crease gradua lly towards higher effective tem- 
peratures, those of iMashonkinal rather display a flat behaviour or 
even decrease. Part of the solution might lie in the adoption of 
different oscillator strengths, placing the lines on slightly different 

parts of the curve-of-growth. 

Of earlier studies we highlight lThevenin & Idiaril jl999h who 
derived metallicities in NLTE for a size-able sample of metal-poor 
and metal-rich dwarfs and subgiants. They report a strong tendency 
of increasing NLTE effects with decreasing metallicity, akin to our 
findings but of greater magnitude. The reason for the difference 
is likely due to different model atoms; their atom lacked e.g. the 
important high excitation levels mentioned above. 

3.5 Conclusions 

The most important NLTE effect for the determination of spectro- 
scopic stellar parameters of late-type stars is the perturbed ionisa- 
tion equilibrium of Fe I and Fe II lines. Departures from NLTE are 
significant for the derivation of Fe abundances from neutral lines, 
at the level of <0. 1 dex in solar-metallicity stars, and < 0.5 dex 
in metal-poor stars, while mostly insignificant for Fe II lines. Sur- 
face gravities obtained through the ionisation balance of neutral and 
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singly ionised lines are consequently affected by 2.5 - 3 times the 
typical NLTE abundance correction, in logarithmic units. Cooler 
stars are less affected than hotter, and dwarfs are less affected than 
giants, in a well-behaved pattern that can be understood by the vari- 
ation of radiative and collisional rates over the HR diagram. 

Due to th e demonstrated impact of <3D> effects on the Fe 
line formation l lBergemann et alj|2012l) . we advice particular cau- 
tion with temperatures derived through excitation equilibrium of 
metal-poor stars. ID LTE effective temperatures are likely to be 
underestimated in this regime and ID NLTE effective temperature 
possibly even more so. However, except for extremely metal-poor 
stars and high-temperature giants, the NLTE effects on the exci- 
tation balance in ID are typically < 50 K. Microturbulence values 
determined from Fe I lines are likewise very little affected by NLTE 
for these stars. 

Our grid calculations can be used to infer individual NLTE 
abundances for thousands of lines in late-type stellar spectra. It may 
be readily predicted that the field of Galactic archaeology subse- 
quently will be confronted with a minor compression of metallicity 
distribution functions, due to the strictly increasing NLTE effects at 
lower metallicity. For example, a distribution function of the Galac- 
tic halo made up of T s g = 5000K/logg = 2.0 giants would be 
affected by +0.3 dex at the hyper metal-poor end, while the metal- 
rich end would remain essentially unchanged. On the other hand, 
the relative abundances of dwarf and giant stars may well be little 
affected, since dwarf samples are commonly hotter than giant sam- 
ples, and consequently they are prone to departures from LTE of 
similar sizes. For example, the NLTE corrections of giants of the 
type just mentioned follow essentially the same behaviour as those 
of T eff = 6500 K/logg = 4.0 dwarfs. 

The propagated NLTE effects on Fel lines on spectroscopic 
parameters are large enough to affect the classification of evolution- 
ary states, ages, and masses. Forced with the requirement to meet 
ionisation balance of Fe I and Fe II, LTE analysis will underesti- 
mate surface gravities and/or overestimate effective temperatures, 
hence moving stars either upwards or leftwards in the HR diagram. 
Finally, trends of element ratios, [X/Fe], with respect to [Fe/H] are 
likely to undergo a change of shape. Assuming no change in el- 
ement X, negative slopes will become flatter and positive slopes 
steeper. 

In the end, accurate Fe abundances is in general a pre-requisite 
for the accurate determination of abundances of all other elements. 
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